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Abstract 



We present a photonic wave packet construction which is immune against the decoherence 
effects induced by the action of the Lorentz group. The amplitudes of a pure quantum state 
representing the wave packet remain invariant irrespective of the reference frame into which the 
wave packet has been transformed. Transmitted information is encoded in the helicity degrees 
of freedom of two correlated momentum modes. The helicity encoding is considered to be par- 
, ticularly suitable for free-space communication. The integral part of the story is information 

■ retrieval on the receiver's side. We employed probably the simplest possible helicity (polariza- 

' tion) projection measurement originally studied by Peres and Terno. Remarkably, the same 

0^ ■ conditions ensuring the invariance of the wave packet also guarantee perfect distinguishability 

in the process of measuring the helicity. 

O 

1 Introduction 

On" 

The helicity density matrix of a generic momentum-helicity wave packet under the action of the 
Lorentz group is ill-defined. The underlying reason is the explicit helicity dependence on the mo- 
mentum eigenstates when transformed from one reference frame into another. Consequently, tracing 
over momentum is an act of violence against the rules of quantum mechanics because one intends 
to sum helicity density matrices from different Hilbert spaces. It makes no sense to talk about 
transformation properties of the helicity density matrices and there is no group representation the 
helicity matrix can transform according to. The situation is further worsened by attempts to mea- 
sure such a state to recover the encoded information. This problem was first systematically studied 
in [I] but some explicit comments were already made in [2] and probably even before. A partial 
remedy was found in pQ too where the authors first realized the importance of a specific helic- 
ity measurement and then introduced the concept of a three-dimensional helicity density matrix 
to deal with the problems mentioned above. Regarding the helicity measurement, it is a simple 
projective measurement onto the plane perpendicular to the direction of propagation of the wave 
packet and due to its simplicity it will be used in this work as well. We will refer to this kind of 
measurement as the Peres- Terno (measurement) scheme. We will, however, abandon the concept of 
three-dimensional helicity density matrix. The additional longitudinal degree of freedom introduced 
in [I] for the sake of having at least definite transformation properties under ordinary rotations will 
be of no use here. Moreover, the introduction of a longitudinal degree of freedom sounds like there 
is something forbidden going on but it is not the case. This expression, however, simply refers to 
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the third degree of freedom of a specific POVM measurement. One could in principle have POVMs 
with arbitrarily many outcomes to be able to reconstruct the wave packet to an arbitrary precision. 
The whole problem of the helicity density matrix has been further clarified and generalized to other 
(non-relativistic) situations in [3] and [5] . Problems of density matrix construction were considered 
in other scenarios, too [6]. For the purpose of this article we stick to the relativistic context. 

There exists a class of wave packets where the definition of the helicity density matrices causes 
no troubles. These are linear polarization wave packets in which case, loosely speaking, all the 
momentum three- vectors point in the same direction (the direction of propagation). This class of 
wave packets has very simple transformation properties and allows for introduction of the invariant 
helicity density matrix [9l [7]. That is, tracing over momentum is allowed and by invariance we 
mean that the coefficients of the resulting density matrix are independent of the reference frame. 
States transforming covariantly under a general action of the Lorentz group are a subject of interest. 
The classification of two-photon states has been provided in Ref. [10] . Let us emphasize that contrary 
to our paper only states with sharp momenta were considered so the physical picture considerably 
differs. 

The aim of this paper is twofold. In the first part we construct localized wave packets where the 
information is encoded in the helicity degrees of freedom (Sec. [5]). The effect of an arbitrary Lorentz 
transformation is studied and the conditions for invariant transformation of these wave packets are 
found. The helicity density matrix derived from these states is well defined and is invariant under the 
general action of the Lorentz group. Finally, the effect of the Peres- Terno helicity measurement on 
the receiver's side is investigated concluding the perfect discriminability of two initially orthogonal 
states in contrast to the original solution. Appendix A contains some further technical details. 

The second part is solely aimed at non-experts on the problem of the photonic wave packets 
construction and their relativistic transformations. It consists of Appendices B and C where we first 
briefly review Wigner's phenomenal contribution to the understanding of the role which is played 
by the Poincare group in quantum field theory in Minkowski spacetime . Appendix C contains a 
discussion of a rather complicated issue of measurement for generic photonic wave packets and the 
solution of this problem found in [T] . 

2 Relativist ically invariant photonic wave packets 
2.1 Photonic wave packets 

We adopt the following notation throughout the article. The components of a general momentum 
four- vector (briefly four- momentum) p is written as p**. We recall that indices from the Greek 
alphabet are reserved for labeling all four components of the four-momentum. The four-momenta 
p 11 living in the four-momentum space can be of null {p^p^ = 0), timelike (p^p^ > 0) or spacelike 
{p^P 11 < 0) character (omitting the zero vector). For the purpose of this article we are interested 
in the first class. Hence we set c = 1 and normalize energy such that p° = 1 so we get pip 1 = 1 
(i € {1,2,3} = {x,y,z}) for all null four-momenta from the future light cone. For further details 
about the Poincare group and how it acts the reader can take a look at Appendix B. Unless explicitly 
mentioned the manifold we will consider in this paper is the four-momentum space and not the usual 
Minkowski spacetime (one temporal and three spatial coordinates). The difference is again sketched 
in Appendix B. For a certainly better exposition, one might also consult extensive literature on the 
subject of the representations of the Poincare group, for example, in [T2] H31 H71 IT5] . 
A generic single-particle pure photonic wave packet can be written in the form 



where a p ,/3 p e C, |a p | 2 + \/3 p \ 2 = 1, / d[i(p)\f(p)\ 2 = l,dju(p) = l/(2(27r) 3 )d 3 j? is the Lorentz invari- 
ant measure and f(p) is an envelope (scalar) function therefore staying invariant under the Lorentz 
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transformation [Tj5]. The helicity eigenstates span a two-dimensional complex Hilbert space and 
these degrees of freedom are perfectly suited for transmission of quantum or classical information 
in free space. We keep the zeroth component of all four-momentum vectors equal to one (or equiva- 
lently that pip 1 = 1 holds) so we restrict ourselves to monochromatic wave packets. We did not make 
our considerations less general by restricting to monochromatic wave packets. The measurement 
introduced later will simply be color-blind, that is, it will act trivially on the color subspace C p for 
every direction p. 

Furthermore, we naturally assume that a, /3 are identical for all momentum directions so we 
may drop the subscript. The preparation of the wave packet is assumed to be under our control. 
Moreover, it would be a strange way of encoding quantum information into the wave packet if the 
quantum state differed for different momentum four-vectors already at the preparation stage. How- 
ever, the momentum directions (the envelope function f{p)) continue to be distributed arbitrarily. 
We thus rewrite Eq. (TIJ as 



U(R(0,<p)) W=a dM(«)/(«)fflcp[*0p,,] \q,+)+0 / dfx(q)f{q) exp [-i9 p , q ] \q,-), (3) 



where i?(i9, tp) determines the spatial rotation of the coordinate system (or of the wave packet if we 
adopt the active point of view) in the direction given by $ and tp and from the scalar character of 
the envelope function it is understood that f(q) = f(p). The transformation in Eq. d3j is indeed 
nothing else than a U(l) transformation and the parameter is the Wigner angle from Eq. (|B-8j) . 
Also, from Eq. (|B-8|) follows that the Wigner angle is explicitly dependent on p and so every ket 
under the above integral acquires a different phase. As pointed out by the authors of p] when we 
trace over the non-interesting degrees of freedom (the momentum) the resulting two-dimensional 
helicity density matrix does not transform according to a representation of the SU(2) group. It does 
not actually have transformation properties at all. We can get an insight into why tracing over the 
momentum is ill-defined. Take just two different directions pi,P2 in Eq. ([2]). The rotation of the 
wave packet given by R("d, (fi) results in an unequal phase change and thus the coefficients originally 
equal to each other will generally become different. So by a suitable rotation two vectors with 
the same components can be made orthogonal and viceversa. Note, however, that even before the 
rotation the operation of tracing over the momentum is not valid despite of bringing up a reasonable 
density matrix. 

2.2 Assembling of invariant photonic wave packets 

We will take a slightly different path towards a reasonable wave packet density matrix. While we 
stick to the helicity measurement studied by Peres and Terno [I] we change the building block from 
which a wave packet is constructed. The choice of the Peres- Terno measurement scheme is well 
motivated. This helicity measurement is probably the simplest one possible. We recall the whole 
issue for readers unfamiliar with the problem in Appendix C. The fixed projection measurement 
breaks the rotational symmetry. This is the reason why we cannot hope in a construction of a 
covariant density matrix. 

Nevertheless, we will be able to construct an invariant density matrix. That is a matrix whose 
components will stay constant independently on the receivers reference frame. How do we achieve 
that? Contrary to [1] we use helicity-momentum entangled two-mode states as the construction 
blocks. We let the helicity-momentum entangled two-mode states propagate in the z direction in 
a fixed coordinate system. We first find the conditions for the entangled states under which they 
do not acquire any phase in a different (spatially rotated) reference frame. Having these conditions 
imposed, we investigate the effect of the helicity measurement on the entangled states and also see 





Let us see what happens if we first just rotate the wave packet 
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the effect of a boost in the z direction. Based on the results, we shall construct wave packets which 
can perfectly transfer quantum or classical information even in the relativistic regime. 




2.3 The effect of spatial rotations 

For the purpose of the helicity wave packet construction introduced later we will now show that 
for the Bell states (Eq. ([7])) with specifically correlated momenta the effect of a spatial rotation is 
harmless. 

Following Appendix B, let us introduce the helicity eigenstates of the standard momentum k 

|-> fc =(lV (4) 

The axial component can be disregard as a result of the application of the Coulomb gauge which 
explicitly sets the component to zero. The helicity eigenstates can be transformed into the linear 
polarization eigenstates corresponding to the standard momentum k by means of the following 
unitary operation 

^ 2 \0 V2j 

We will need this operation to transform the rotation operator from the linear polarization basis 
R(d,ip) [T] to the helicity basis R(d,ip). Therefore, to determine the helicity vector for an arbi- 
trary direction parametrized by two angles ip we transform the helicity eigenstates as |±) p = 

R(d,v)\±) k = SR(d,v)S- 1 \±) k . Then 

^ /(cos ■& + 1) exp (— iip) (cos $ — 1) exp (— iip) \/2sin-# exp (— iip)\ 
-R(#, ip) = — I (cosz? — 1) exp (iip) (cosfl + 1) exp (iip) v^sin^exp (iip) (6) 
V -v^sintf -V2sintf 2cosi? / 

and so throughout the article we work strictly in the helicity basis 0. 
We write down the usual four Bell states 

|*(P1,P2)> ± = ^= (bi, +) \P2, +) ± \pi, -) \P2, -)) (7a) 

|*(pi,pa)) ± = ^=(|pi,+)b 2! -}±bi ! -)b2,+», (7b) 

where \pt , ±) = (pi, ±) \vac) using the notation of Appendix B (from now on, the Latin subscripts 
label corresponding four- vectors and not their components) . No need to repeat that the states above 
are not wave packets yet. Using Eqs. ([5]) and ^ we rewrite states (|7a|) and (|7b|) . The helicity triads 



1 Note that we get the helicity eigenstates in the linear polarization basis by the action of S 1 |±} fc . Before 
proceeding, let us mention that there exist many conventions for the matrix S. One can, for example, often find 

/-ll 
S' = — -i -i 



v2 \ \/2, 

The rotation matrix then looks different as well as other objects we will meet later but the physical consequences are 
equivalent. By working with Eq. J5} we employed the convention of 1121 . 
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then read 



(cos$ + 1) exp (—i<p)^ 
\P,+) = $ ( (cos 0-1) exp (ty>) j (8) 
-V^sintf 



1 



'(cos $ — 1) exp (—iip)^ 



\P>~) = o (cos 0+1) exp (ip) . (9) 
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An arbitrary rotation can be decomposed into rotations around the z axis followed by another 
rotation around the y axis. Following the prescription sketched in Remark Q] (or looking into [IB] 
or eventually [9] where the Wigner angle was calculated in the representation of the cover of the 
Lorentz group - the group SL(2, C)) we find that the rotation i? z (A) around the z axis through angle 
A induces the Wigner phase angle 9 z ^ p = 9 p = for all p as long as p ^ k. Exactly the measure zero 
possibility p — k is excluded from the definition of \<S>(jp\,p2)) because that would imply p\ = pi = k 
and |$(pi,P2)) to be converted into \<j>) oc ([o^fc, +)] 2 + [at(k, — )] 2 ) \vac) due to indistinguishability 
of such photons. Hence the Bell states transform as 

U{R Z (\)) : |*(pi,P2)) ± -► |$(i?,Pi,i? z p 2 )> ± (10) 

U(R Z (X)) : |*(pi,p 2 ))± -► |*(i?,Pi,i? z p 2 )) ± • (11) 

The rotation R y (w) around the y axis results in a more complicated formula for the Wigner angle 

" sin rosing 1 

sin w cos $ cos ty? + cos w sin ?9 

To tackle the momentum dependence we will assume the momenta p\ , pi to be correlated such that 
the net effect of the rotation will be the overall phase equal to zero for all p\. That implies for 
Eq. (|7a|) to find such pi that the corresponding Wigner phases are conjugated (because the helicity 
is equal for each product). Having = (1, sin)?; cos^j, sini^ sin ipi, cosi?i) this condition implies 
that we have to find solutions of the equation 

sintnsin^i sinn7sin(<^i + x) 



9 p = arctan 



sin vj cos $i cos ip\ + cos w sin -Q\ sin w cos($i + y) cos(<y9i + x) + cos w sin(i?i + y) ' 



(13) 



where <pi + x — (f2, i?i + y = This equation deserves an explanation. To find a solution means 
to find x, y independent on m. In other words, we want to locate all fixed points. Note that these 
points are not completely fixed since they vary with ^i,i9j.. 

Similarly, for the second couple Eq. (|7bp we have to find pi , pi such that the Wigner phases are 
equal (the helicity is opposite for each product) 

sinccjsintp! sin w sin(<^i + x) 



sin w cos $i cos (p± + cos w sin d\ sin w cos(z?i + y) cos{ip\ + x) + cos w sin(i?i + y) 



Lemma 1. There exist four fixed points for Eqs. M3\) and Jj^| ). 
Proof. Eq. ([IB")) can be cast into the form 

tan zu [snupi cos (d\ + y) cos (epi + x) + sin (ipi + x) cos i?i cos ifi] = — sin (ipi + x) sin#i — sin^i sin ($i + y). 

(15) 

To have solutions independent on w and satisfying the above equation the only possibility is when 
both sides are equal to zero. That implies 

cos + y) cos ((pi + x) cos i?i cost^i 

: — 7 \ = : (16) 

sin (ipi + x) sin i^i 

sin (gi + y) _ sin^i 

sin (ipi + x) sin ipi 



5 




Figure 1: Four fixed points corresponding to Figure 2: Here we demonstrate how the po- 

Eqs. (|18aj) . (|18bl) . (|18cp and (| 18d[) are indicated sition of fixed points changes for a different 

by points a, b, c and d. Curves corresponding to momentum direction. The fixed points sat- 

a different rotation angle w about the y axis will isfy Eq. (fT4f so in contrast to Fig. Q] there is a 

always meet in these points. The fixed points sat- global (trivial) fixed point in the centre which is 

isfy Eq. (fT3|) . One can notice several gaps on every constant for all 0i, ip\. One can notice several 

curve. These are singularities of Eq. (fl3|) . gaps on every curve. These are singularities of 

Eq. (HI. 



It is easy to enumerate all possibilities for which, for example, the first equation is satisfied and then 
to check if it also holds for the second one. We end up with four different solutions (modulo 2tv) 

<$i = —<pi, $2 = 0i (18a) 
ip 2 = tt - <pi,i?2 = — #1 (18b) 
<P2=ir + tpi,'&2=w — l 8i (18c) 

<p 2 = <Pl,$2 = 01 +7T. (18d) 



Having the solutions of Eq. (|13|) we immediately get the solutions for Eq. (|T4)) . The later equation 
is obtained from the former one by a mere shift x i— > x, y i— > y + ir. Hence 

V>2 = -^i,i?2 = +7T (19a) 
^ 2 = 7r - <^i,02 = tt - 01 (19b) 

!^2 = 7T + </5i,02 = — 01 (19c) 
(f2 = <£l, 02 = 01- (19d) 

In Figs. [I] and [5] we see the position of four fixed points for some 0i, (fx. If we had in Fig. [5] the 
same parameters 0i,<^i as in Fig. [1] the picture would indeed be only shifted along the y axis. We 
can check that all the solutions satisfy the corresponding Eqs. (fT5)l and (fH|) . ■ 

The first solution Eq. (|I 8aj) is the most interesting one since p2 dwells in the same hemisphere as 
p\. Commenting on the rest of the solutions, Eq. (|18b[) is equivalent to Eq. (If 8ap and the remaining 
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two solutions are a bit awkward because p 2 points in the opposite direction. We might, however, 
say that 

U{R v {w)):\$(p l ,p 2 )) ± ^\${R y p 1 ,R y p 2 )) ± . (20) 

From the physical point of view there is a noteworthy asymmetry in the solutions for Eq. (|14p. 
The first two fixed points again correspond to a wave packet with the second momentum pointing in 
the opposite direction and the last two solutions are just an identity operation (p 2 = p\ ) . So there 
are two solutions where 

U(R y (w)) : |*(pi,P2)) ± -> \*(R y p 1 ,R y p 2 )) ± (21) 

holds but, for the practical purposes, the opposite momentum direction is an obvious obstacle. 

The conclusion is that as long as we keep the momenta in the Bell states correlated as described 
above it is not necessary to introduce three-dimensional helicity density matrices since there is no 
phase change due to the spatial rotation at all (cf. Eqs. (ITU1) . (|TTj) . (|2U)l and (|2"Tj) ). 

The way of getting rid of the Wigner phase for spatial rotations is reminiscent of the method 
presented in [19j . The authors constructed wave packets by entangling two spatially distinguishable 
wave packets carrying the same momentum and opposite helicities. Every spatial rotation induces 
an opposite Wigner phase for each wave packet so the resulting phase is zero 0. As we already 
stressed, however, the Bell states we work with so far are not wave packets so the construction 
presented here is qualitatively different. 

2.4 The effect of Peres- Terno measurement 

The projection present in the Peres- Terno measurement process (Eq. (|C-2jl ) acts by 'cutting off' the 
z component of Eqs. (jHJ) and ([9]). Henceforth, after the proper normalization we get 

IK • Id +)-> I v +) ~ - ( (C ° S d + 1} CXP (~ i<P A (22) 
n fe . |p, +; -> [p, +) - N y (cosi9 _ x) exp ^ j W 

life ■ \p, ) -> IP, )- N y ( cos tf + i) exp (itp) J ' [Z6 > 

where N = v2 + 2 cos 2 Note that (p, +\p, — ) = for all p but (p, +\_p, — ) ^ as expected. The 
'floored' kets indicate the action of the projection 11^ on a momentum-helicity state. 

Interesting things start to happen when we ask how the projection Ilfc acts on the Bell states. 
First of all, these states are two-mode states so the projection is actually III 1 ' ® IT^ 2 '' where the 
superscripts indicate which mode is being measured. The action of the projection can be found in 
Appendix A. In general, these states are subnormalized and non-orthogonal and thus resembling a 
general situation studied in pQ. But, intriguingly, when plugging in some of the solutions for the 
momentum correlations (namely Eqs. (|18aj) . (|18c[) . (|18dj) . (|19aj) and (|19bjl ) found from the study 
of fixed points for arbitrary rotations about the y axis we find that (i) the resulting floored Bell 
states are orthogonal for all and (ii) one of the states is always (up to a normalization) 

invariant. This has some interesting consequences particularly for the solution Eq. (|18al) . Recall 
that this is the only non-trivial solution where the resulting state 'points' in the same direction 
as the propagation direction. From now on we will stick just to this solution since it will be later 
relevant for our discussion of localized wave packets (the analysis for the rest of the non-trivial 
solutions gives similar results but their use as localized wave packets is obviously none). 

When the previous claim written in detail for ip2 = —<pi,'&2 = $i we get from Eqs. (|A-2|) 

2 Note that this construction is not working for arbitrary wave packets. It works only for wave packets where the 
rotation induces the same Wigner phase for all momenta from which the wave packet is 'assembled'. 
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and fO|) 



n 



(i) 



n£> : |*(p 1>P2 )>. 



L*(pi,p 2 )>+ = 



i \ 

exp (2^)^1=1 

1 / 



normalization 



> ai L*+) +a 2 L*+) +™ 3 L*-) (24) 



n 



(i) 



ni 2) : |*(pi,Pa)>_ -> L*(Pi>P2))_ = 



1 



/ 2cosi?i \ 



\-2cos??i/ 
l/V2([0i) L0 2 > + Ui 



1 + COS 2 ??1 



normalization 



(25) 



LI2)) and similarly for the 



where m G R and £V a? = 1. We denoted 
rest of the states. Note that both output states are orthogonal for all i?i and <p\. Also, the resulting 
states are subnormalized even though states (f22|) and ((23)) are normalized. This is again the effect 
of cutting off (alias projection) the 'longitudinal' part of Eqs. (jHJ and (J9j> . The normalization factor 
for Eq. (El and Eq. 051) is - i* 0082 f 1 ,, and 1 + cos2 f 1 , respectively. Recall that the measurement 

. We also want to draw attention 



2\/l+cos 4 tfi 

indeed acts as a trace-decreasing non-completely positive map 



to the similarity between the behavior of the investigated Bell states under the irrep of the product 



SU{2) x SU(2) and the projection IT^ f^n^. In the fist case the singlet is preserved and the triplet 
states transform among themselves. This is a starting point to considerations about reference frames 
and decoherence-free subspaces |20j . In our case one of the triplets is invariant meanwhile the second 
triplet mingles with the singlet and the remaining triplet state. 

An important comment also related to decoherence-free subspaces has to made. As we have 
seen, the Bell states remain orthogonal after the projection 



T (2) 



(*(P1,J)2) + |*(P1,P2)-) =0 



($(pi,P2)+l*(pi,P2)-> = 0. 



However, the action of the projection is not SU(2) x SU(2) covariant. In other words, there is a 
preferred basis spanned by states from Eq. ([7]). One faces a slightly analogous situation encountered 
in the study of decoherence-free subspaces [21]. However, the presented results still hold for any 
orthogonal basis \pi, +), \pi, — ) forming the actual entangled states in Eq. ([7]). The only difference 
is given by the fact that our underlying Hilbert space is infinite-dimensional. 



2.5 Boost along the z axis followed by measurement 

The found result has several serious consequences. First, if we prepare two orthogonal pure qubits 
having the basis formed by the Bell pair |$(pi,f>2))±, they will be perfectly distinguishable by the 
measurement. This was not the case for the original Peres- Terno scheme. But more important, if 
we boost along the z axis and then perform the same projective measurement we see that state (1241) 
remains invariant in the sense that the boost only modifies i?i(= #2) and keeps ipi(= —^2) intact. 
Of course, by the change of $2 the four- momenta pi,P2 change too as Eq. (|B-1[) dictates. They 
symmetrically close in or open up depending on the observer's velocity direction. In a real experiment 
this would be registered as the color change (this is the Doppler effect which we formally suppress 
by the renormalization of the new four- vectors) . 

Hence the effect of a boost in the z direction B z (ri) followed by the projective measurement 

<g> n^ 2) is 

nj> x) ® Hf } : \$(B zPl , B z p 2 )) + -> L*(-B,pi, B zP2 )) + (26) 
IL{ 1] ® nf > : \$(B z pi,B zP2 ))_ -> mB zPl ,B zP2 ))_ . (27) 
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Let us conclude this section with two comments, (i) As in the previous case all output states are 
orthogonal for all boosts and every p±. We have to stress, however, that the above transformations 
are not unitary. We recall that the normalization of Eqs. (l22l) - (|25p is explicitly i?— dependent. It 
would be a clear sign of inconsistency if we had a finite-dimensional unitary representation for boost 
transformations - the boost generators are the noncompact part of the Lorentz algebra. Similar 
inconveniences were encountered in pQ but in our case it does not cause serious troubles, (ii) For a 
state propagating in the z direction the most general Lorentz transformation can be expressed as 
R z (X)Ry(zu)B z (r]). Hence, we have shown that the Bell states propagating in the z direction and 
with the appropriately correlated momenta are unaffected by an arbitrary rotation (and therefore 
a general boost) if measured according to the Peres- Terno scheme. 

2.6 Relativistically invariant photonic wave packets 

The real touchstone for quality will be the behavior of wave packets composed from the previously 
studied entangled states under the Lorentz transformation. As indicated, we will demonstrate the 
construction and properties of invariant wave packets with \&(pi,P2)}± serving as the logical basis 
for construction of a (pure) wave packet qubit. We will not have anything specific to say about the 
actual feasibility of producing of such states in a laboratory. Also, the dispersion effects of other 
than relativistic origins are not considered here. 

Let us follow the prescription and notation for the single-mode case Eq. (fT]) 

|0 12 ) = J d f i(p 1 )dfi( P2 )f{p 1 , P2 )(a Pu2 |*(pi,p2)> + |*(pi,fti)>_). (28) 

Similarly to the situation of ordinary wave packets we assume a and (3 to be momentum independent 
(see the discussion below Eq. (j2J). 

Taking into account Eqs. (|10[) and ([20)) and how the boost acts we easily find 

|Ui2> = U{R z {\)R y {m)B z (r,)) ^12) 

= a / d/i(gi)dy*(ft8)/(?i,«a)|$(gi,ft)) + +j8 / d/i(«i)d/i(«2)/(gi, q 2 ) &))_, (29) 



where qt = R z {\)R y {m)B z (r))pi. As discussed earlier in Eq. ([3]) the measure and the envelope 
function are Lorentz invariant (/((ft, 52) = f(pi,P2))- Hence, we call the transformed wave packet 
invariant because the state amplitudes remain unaltered by an arbitrary Lorentz transformation. In 
other words, otherwise disastrous decoherence between helicity and momentum degrees of freedom 
is kept under control by a careful choice of the logical basis. 

As a final step, the measurement described by Ilj^ <8>nj^ is realized, due to Eqs. (f24| and (|25|) . 
by the projections 

r 1 = L$-X*-J, r 2 = i-r x . (30) 

To be more precise, we have seen that the map related to the measurement operator n[ (g) n£ is 
trace-decreasing. Consequently, this measurement is a conditional (or post-selected) measurement. 
The overall probability of measurement is lower than one and depends on the details of the wave 
packet, namely on the choice of its envelope function. When the measurement event occurs the 
observer can always perfectly distinguish both orthogonal basis states since the two subspaces from 
Eq. (f30]) are orthogonal and independent on the observer's reference frame. Note that we tacitly 
but naturally assume the existence of a global coordinate system. 

3 Conclusions 

There exists a class of photonic states we call the realistic wave packets. It is a weighted and nor- 
malized superposition of common eigenvectors of the momentum and helicity operators. The reason 
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we call them realistic is twofold. From the mathematical point of view, if the corresponding enve- 
lope function is square integrable (but otherwise almost arbitrary) function, such an object lives in a 
Hilbert space. From the physical point of view, this object can in principle be prepared, transmitted, 
manipulated in a localized manner and finally detected in a laboratory. The problem appears when 
we consider wave packets in a relativistic regime. The general Lorentz transformation related to the 
change of a reference frame acts individually on every momentum and induces momentum-helicity 
entanglement. The situation is even worse when we try to recover the information by an act of 
measurement. The simplest (and probably the only feasible) way is the Peres- Terno scheme of the 
helicity projection onto the plane perpendicular to the direction of propagation. 

Surprisingly, in this paper we have shown that under these unfavorable conditions we are still able 
to prepare localized wave packets for which the information can be, at least in principle, transmitted 
and recovered in every reference frame with perfect fidelity. There are two key ingredients. First, 
the logical qubit basis are two-mode maximally entangled helicity states. Second, the momentum 
degrees of freedom of these states are correlated in a precise manner. This results in canceling 
the explicit momentum dependence for an arbitrary rotation and boost. Intriguingly, the same 
condition implies that during the measurement according to the above scheme the orthogonal logical 
basis states forming a qubit are projected into two orthogonal subspaces enabling us to perfectly 
distinguish them. 

It has been known from the previous analysis that the map governing the evolution and mea- 
surement of generic wave packets in the relativistic setting is trace-decreasing and non-complctely 
positive. It is a consequence of the non-covariant helicity measurement scheme and this leads to 
all aforementioned effects and troubles. Here we have encountered precisely the same behavior 
but contrary to the previous work these effects are now harmless. Consequently, using the investi- 
gated encoding the problem of non-covariant transformation of the von Neumann entropy of helicity 
density matrices disappears. 

Finally, let us stress that the choice of the Peres- Terno measurement (arguably spoiling all nice 
transformation properties of the momentum-helicity eigenstates) is directly motivated by its relative 
feasibility in a laboratory. Should more sophisticated helicity measurements be available we might 
obtain more elegant results. 
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Appendix A 

The projection corresponding to the Peres- Terno measurement scheme acts as 



(A-l) 



where the 'floored' quantities indicate the projection action. Denoting 



1 



1 



N V / (l + COS 2 l?i)(l + COS 2 02 ) 



we get 




$(pi,P2))+ L$(pi,p 2 )) + 



N 



1 



/ (costfi cos $2 + 1) exp (—i((pi + <P2))\ 
(cosi?i cosi9 2 — 1) exp (— i(ip\ — <f 2 )) 
(cosi?i cos $2 — 1) exp (i(<pi — ^2)) 
y (cos^i cos??2 + 1) exp (i(ipi + ^2)) / 



(A-2) 




N 



1 



/ (cos^i + cosi9 2 ) exp (—i(<pi + </>2))\ 
— (cosi?i — cosi?2) exp (i((f2 — fi)) 
(cosz?i — cos$2) exp (—i(ip2 — <fi j) 

\-(cos$i + cosz? 2 ) exp (i(<pi + ip 2 )) / 



(A-3) 




/ (cos^i cosi? 2 - 1) exp {-i(tpi + </?2))\ 
(cos^i cosi? 2 + 1) exp (i(ip2 - 
(cosi?i cosi? 2 + 1) exp (—i(if2 — (pi)) 
\ (cos^i cosi? 2 - 1) exp (i(<pi + 1^2)) J 



(A-4) 
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and 



n« ® 4 2) : |*0i, P2 )>_ -> L*(pi,P2)) 

Appendix B - Review the unitary representation of the Poincare 
group 

The purpose of the following two Appendices is to recall the problem of the polarization encoding 
for realistic photonic wave packets. Readers not familiar with this problem full of intricate details 
might find this short review useful. 

The Poincare algebra has two Casimir operators C\ — P X P X and C2 — W X W X where W x = 
— 1/2 e xpau P p J au] is the Pauli-Lubanski vector, { J23, J31, J12} = {Ji,J2,Js} are the total angular 
momentum generators and { Jq%, J02, Jm} = {Ki,K2,K 3 } are the boost generators. We recall 
that e xpau ' is the Levi-Civita tensor and e 0123 = — 1. The corresponding Poincare group induces 
translations x^ — > h^x v + (in the coordinate representation) where A(£) £ SO(3, 1) is a Lorentz 
transformation A(£) = exp (— i/2C <J? Jo-?) and T(a) = exp (— ia^P^). It is well known that the 
Poincare group is a semidirect product of the proper orthochronous Lorentz group and a group of 
translations SO(3, 1) k R 4 . 

The group acts on set K which is a linear vector space equipped with the Minkowski metric with 
the signature { — h ++}■ Hence the set K is a manifold rich in structure. Note that in the case of 
Minkowski spacetime the action of the Poincare group is uninteresting since it acts transitively. It is 
not so in the four-momentum space. The Lorentz group acts in a similar manner for both spaces but 
the action of the translation operator differs. One can see it easily (in an admittedly handwaving 
way) realizing that the Fourier transform of a function / : Xi 1— > xi — a% for some constant a is a 
simple phase transformation in the three-momentum space. 

The analysis of the four-momentum space K shows that there are six orbits. The whole classi- 
fication from the physical point of view can be found, for instance, in |12j 113] . We are interested 
in the orbit which physically corresponds to null vectors with a positive zero component of the 
momentum four-vector. They correspond (after quantization) to the scarce but important class of 
free massless particles containing photons and possibly gravitons. By the very definition of orbit no 
Lorentz transformation can change the character of a null vector. In other words, photons always 
travel at the speed of light no matter the reference frame a potential observer resides in. We are free 
to generate the 'null orbit' from any null momentum four-vector but it is customary and advanta- 
geous to choose the simplest possible one - the standard four-momentum — (1, 0, 0, 1). Note that 
other choices are completely equivalent as discussed in detail in [17]. From the physical perspective, 
the standard direction corresponds to (eventually quantized) plane electromagnetic waves (photons) 
traveling at the normalized speed of light along the z axis. Nevertheless, the word 'physical' should 
be used carefully. It is not, in principle, possible to generate such a field (classical or quantum). 
The reason is that the sharp momentum value implies that the classical wave is 'infinitely' spread in 
position and one would need an infinite amount of energy to prepare such a field. In the quantum 
case the corresponding state does not even occupy a Hilbcrt space of square integrable functions 
and is completely delocalized in Minkowski spacetime. Classical or quantum fields created in a 
laboratory are in reality wave packets with a finite spread in momentum and position. 

Using another property of the algebra [P^, P v \ = and C\ = C2 = which holds in the massless 
case only we denote a 'single-particle' state \p, a) as an eigenstate of P M where p p are eigenvalues of 
P**. Other possible degrees of freedom related to rather than to C2 (note that \W P ",P 1J ] = 0) 
are gathered in a 

P»\p,a)=p»\p,a). (B-l) 



1 

N 



/-(cos$i - cos-i? 2 )exp (-i(ipi + (f2))\ 
(cos$i + cos $2) exp (~i(ipi — ^2)) 
-(costfi + cosi9 2 ) exp (i(tpi - ^2)) 
\ (cosz?i - cos-# 2 ) exp (i(<pi + tp 2 )) I 



(A-5) 
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To be able to use quantum mechanics in the relativistic context as we habitually do in ordinary 
quantum mechanics we need to introduce a Hilbert space and a unitary representation in it which 
respects the composition law for the Poincare group 

U({A, b})U({A, a}) = exp (iT)U({AA, Aa + b}). (B-2) 

The appearance of the phase is due to the projective nature of quantum states (rays instead of 
vectors). For simply-connected groups it is possible to get rid of the phase factor (T = 0). This is 
not the case of the Lorentz (or Poincare) group and it will have some impact at a later stage. 

The Poincare group is non-compact and therefore there is no faithful finite-dimensional unitary 
representation. This might have been a huge obstacle for studying the evolution of relativistic 
quantum states but Wigner realized how to circumvent this problem. He found that the unitary 
action of the Lorentz group is governed by the irreps of its stabilizer subgroups 13 . Loosely 
speaking, the little group 'induces' its action on the Lorentz group. The method is called the 
method of induced representations and was later studied and generalized by [2] and others. The 
Wigner method prescribes that the unitary action of a Lorentz transformation A^p 11 = p lv reads 

[/(A) \p, <t}=J2 Da,„> (S k ) \ P ', </) , (B-3) 

a' 

where D aylJ i(Sk) is the irrep of the stabilizer group Sk keeping the standard null direction k intact. 
The stabilizer group can be written 

S k = L- 1 (A P )AL( P ), (B-4) 

where L(p) : k — >• p takes the standard four- vector to any other null vector (i.e., to a vector on the 
same orbit). It turns out that the stabilizer group for the standard null orbit is the Euclidean group 
ISO(2) (also known as E{2)) transforming a two-dimensional (Euclidean) plane into itself. It is a 
three parametric group generated by translations in two directions and a rotation around the axis 
perpendicular to the plane. The components of the Pauli-Lubanski vector are the generators of the 
corresponding Lie algebra. The group is noncompact so it might seem that we did not improve our 
situation but we actually did. To make the long story short [12] it appears that the noncompact part 
of the stabilizer group is not physically relevant to the evolution of the photon under the Lorentz 
group so the only SO(2) subgroup of ISO(2) remains. The corresponding algebra generator is 
W° = W 3 = J 3 and we find 

J 3 \k,a) =a\k,a). (B-5) 

We can see from Eq. (|B-5j) that 

U(R(ti, cp))J 3 U^(R(d, <p)) \p, a)=a \p, a) (B-6) 

so we define hclicity as the projection of the angular momentum along the direction of motion and 
a is clearly a relativistic invariant. 

The SO{2) group element is as usual recovered by exponentiation R z (9) = exp (— i9 J 3 ) so for 
the purpose of Eq. (|B-3[) we get 

D a , a > (R z (6)) = exp [-ia0]5^ (B-7) 

and so for an arbitrary Lorentz group A : p — > q 

17(A) \p, a) = exp [-ia6 p J \q, a) . (B-8) 

Remark 1. We will call 9 p . q the Wigner angle and it can be explicitly calculated for a given Lorentz 
transformation A : p — > q from Eq. (_B-^[ ) by setting R z (6 Pt q) — Sk- Notice that the phase factor 
exp [— ia9 Ptq ] g U(l) is indeed a unitary representation of R z (9 Pi q) G SO{2). The two groups are 
known to be isomorphic. 
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Remark 2. The parameter a can take just integer or half-integer values as the result of the inability 
to set the phase in Eq. HB-2\) to zero \12f . Moreover, the momentum projection in two opposite 
directions p and —p are related by the parity operator so for photons as massless vector bosons we 
get a = ±1. We see that it is necessary to enlarge the Poincare group by including the parity 
transformation [2]. Clearly, there is no helicity equal to zero since there is no rest frame for a 
photon. 

Remark 3. A state \p, a) coincides with the action of creation operators a'(p, a) \ vac) to the vacuum 
state from the usual quantization procedure of a free massless vector field yielding a single particle 
state of a momentum p and helicity a — ±1 (alias circular polarization). It is straightforward to 
generate multi-particle states by a repeated application of the creation operators for different modes 
satisfying [a(p', ct'), (v{p, <r)] = 8{p' — p)S a ' (T and thus to create the familiar Fock space as a direct 
sum of the completely symmetric Hilbert space of n photons J- = ©.^L TL^ m . 

Also notice that one of the consequences of the Wigner procedure is a clear indication that there 
are just two 'spin' degrees of freedom and they are always perpendicular to the direction of motion 
and to each other. We know this fact without mentioning the Coulomb gauge whatsoever. The 
Coulomb gauge can achieve the same goal but one has to pay the price of not having a manifestly 
covariant gauge condition. That is, the gauge condition must be imposed for every reference frame 
separately. On the other hand, the Coulomb gauge is important and useful in quantum field theory 
from a broader point of view. 

Appendix C - Measurement troubles and wave packet helicity 
density matrices 

The authors of [T] proposed a very simple measurement which corresponds to what actually might 
happen in a laboratory when measuring the helicity (polarization) degrees of freedom. Assume that 
the direction of propagation of the wave packet is the standard direction k. It is customary to place 
a polarization analyzer perpendicular to the direction of motion. Our helicity eigenstates \p, ±) can 
be represented as two orthogonal complex four- vectors £p±- For a free field the Coulomb gauge 
Pi£p ± = implies the axial component of the helicity vector to be zero e° ± = (holds for all p 
because just rotations are considered now). Quantum electrodynamics indicates what happens next. 
We summon up the well-known relation for the helicity eigenvectors |12) 



where the overbar indicates complex conjugation. Following the notation of Eqs. (JSJ) and © and 
Appendix B the helicity measurement for p = k corresponds to the trivial projection 



where special symbols have been assigned for the helicity eigenstates because we work with them 
in the main body. The problem is that the projection is of rank two. When a single photon arrives 
in a direction p ^ k it has a longitudinal component and hence it gets cut during the course of 
measurement. Then, to have a consistent definition of the helicity density matrix in a rotated 
reference frame, one has to reconstruct the longitudinal part because that is the place where some 
parts of the wave function 'got lost'. This basically means to measure the helicity in the planes 
perpendicular to the x and y axes in our fixed coordinate system to be able to build a three- 
dimensional helicity density matrix. Finally, we get a two-dimensional effective density matrix g e f{ 
by cutting a 2 x 2 block from the whole three-dimensional matrix. The resulting state is positive- 
semidefinite and subnormalized. It is the consequence of the lemma following the definition. 




(C-l) 




(C-2) 
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Definition 1. For a matrix of dimension n the leading principal submatrices are all the upper left 
submatrices of dimension k < n. 

Lemma 2. Let g be a density matrix written in terms of (not necessarily orthogonal) pure states \<pi) 
such that g = ^2™—i Pi\<f>i)(4>i\ ■ Then every leading principal submatrix of g is positive semidefinite. 

Proof. Taking the k— dimensional leading principal submatrix of \<pi) means to replace the last n — k 
entries of \<fii) by zeros. We will label the resulting vectors as The corresponding matrices 

\4>i\4>i\ are still positive semidefinite and the positivity stays preserved by taking their convex 
combination g = X)i=i Pi L^«X0iJ ■ The resulting density matrix is subnormalized. ■ 

Corollary. The projection 11^ does exactly the job of cutting out the principal two-dimensional 
submatrix from the 3x3 helicity matrix (recall that the axial part is set to zero for all p by the gauge 
condition). Therefore 



where N is a normalization constant. 

Eq. (1C-3|) says nothing other than if we measure the helicity perpendicular to the propagation 
direction (the z axis in this case and so life is of the simple form Eq. (IC-2[) ) a statistical mixture of 
all x, y helicity components for all p is generated. This is exactly the result from [Tj. Of course, we 
may decide to measure along a general axis g ^ k. In this case Eq. (|C-3[) is still valid if we take 
Ilfc — > Ilg and render the helicity vectors in the new basis. 

Remark 4. One can adopt the approach of [3] where the above Peres-Terno measurement scheme 
was reformulated in terms of measurement of the Stokes parameters of the incoming wave packet. 
We will not go into detail but we just point out the Lie algebraic aspect of this approach. In a 
similar manner as above the authors basically constructed three different effective helicity matrices 
which correspond to the above mentioned measurement of the helicity in three orthogonal spatial 
directions. From these they constructed a three-dimensional helicity density matrix which in general 
can be written in terms of the su(3) Lie algebra generators A| J18\j 



where are the actual measured Stokes parameters. We now realize that the su(3) Lie algebra is 

composed of three mutually dependent su(2) Lie algebras and they precisely form the three effective 
helicity density matrices calculated in 01/ (and one of them is Eq. fC-3\) ). 

Remark 5. We should stress that in the above process one indeed gets a longitudinal part of the 
wave function but there is no need to call it a non-physical situation fl]/. We simply get a helicity 
component in the plane not parallel to the plane where the helicity is measured. In principle, this 
situation cannot be avoided unless one works with single-particle states or very narrow photonic 
wave packets. 

Also, the three-dimensional helicity density matrix concept has only little to do with the bosonic 
nature of the photon. It is just a coincidence that we measure along three perpendicular directions 
corresponding to our usual spatial dimensions. If we had had a more sophisticated measurement 
device we could have measured the helicity in each and every plane perpendicular to every p direction 
and map it to some other (non-photonic) multi-level quantum system. We would avoid projecting 
it down to our usual three-dimensional space where helicity measurement devices usually operate. 
Pushing it to the limit, this way we would have been able to recover the whole density matrix because 
(p, ±|p',±) = 5{p — p') holds in the infinite- dimensional helicity-momentum Hilbert space. 



Qes = -T7 / dn( P )\f( P )\ 2 n k (ap|+L+A>l-)J(fip<+L+j8 P HJ n ^ 




(C-3) 




(C-4) 
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So far we talked just how a wave packet transforms under an ordinary rotation. The remedy 
was to introduce a three-dimensional helicity density matrix or to eventually create a narrow wave 
packet where the 'longitudinal' component is negligible pQ. The effect of a boost on a wave packet 
is more disastrous. Assume again that the wave packet Eq. (jTJ) is propagating along the z axis 
with the standard momentum k and an observer is also moving along the z axis with a constant 
velocity v = v z e (—1, 1). Having rj = atanh v z the boost B z : p — > q where B z {rf) = cxp (—irjK^) 
does not induce any phase change [TB] so the boost transforms a general single-particle state as 
U (B z (rj)) : \p, a) — »• \B z p, a). The effect of the boost on a four- momentum vector is twofold: (i) The 
magnitude of the three-vector and the zeroth component are no longer equal to one. This is the 
Doppler effect we suppress by renormalization as discussed below Eq. (JTJ . (ii) The gauge condition 
is violated since it is not a Lorentz covariant condition and must be imposed for every q as is usually 
done in quantum electrodynamics [12]. Thus B z : e Pj(T —t e q>a is followed by 

Coulomb gauge : e qi „ -t e q . a - g q , a q, (C-5) 

where g q ^ = e qa /q° (q° = 1 following from (i)). The gauge condition Eq. (|C-5|) assures that the 
new helicity four-vectors stay orthogonal in a plane orthogonal to q. This is a final blow to our effort 
to have nice transformation properties of the helicity density matrix between two different frames 
of reference. 

Remark 6. The things get more technically complicated if we allow for a boost in a general direc- 
tion. Then, a Wigner phase appears what can be easily seen if we realize that such a boost can be 
decomposed in a boost in the z direction followed by a rotation in the required direction. 



16 



